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Abstract 

A new measurement of the cosmic ray antiproton-to-proton flux ratio between 1 and 100 GeV is 
presented. The results were obtained with the PAMELA experiment, which was launched into low- 
earth orbit on-board the Rcsurs-DKl satellite on June 15*'^ 2006. During 500 days of data collection 
a total of about 1000 antiprotons have been identified, including 100 above an energy of 20 GeV. 
The high-energy results are a ten-fold improvement in statistics with respect to all previously 
published data. The data follow the trend expected from secondary production calculations and 
significantly constrain contributions from exotic sources, e.g. dark matter particle annihilations. 

PACS numbers: 96.50.sb, 95.35.+d, 95.55.Vj 



2 



Antiprotons can be produced from collisions of energetic cosmic ray particles, primarily 
protons, with the constituents of the interstellar gas such as hydrogen and helium. Possible 
primary sources of galactic antiprotons include the annihilation of dark matter particles 

m 

and the evaporation of primordial black holes 4]. Cosmic ray antiproton experiments can 
probe production and transport properties of cosmic rays in the galaxy and search for ev- 
idence of exotic production mechanisms. However, such detailed studies of the antiproton 
energy spectrum require measurements with good statistics over a large energy range. Cos- 
mic ray antiprotons were first observed in pioneering experiments in the 1970s by Bogomolov 



e. al. a a„d OoMe. et a., u.n, bal.oon-bo.e .a,„e«c .pect.^ete., Bogomolov et a.. 

observed 2 antiprotons in the kinetic energy range 2-5 GeV while Golden et al. observed 28 
antiprotons in the range 5-12 GeV. Several other experiments followed, covering the kinetic 
energy range 0.2-50 GeV. More than 1000 antiprotons have been observed in the kinetic 
energy range 0.2-4 GeV by the BESS expeHjnent Q while the^ati.tic= at Mghet energies i. 
very limited. The CAPR1CE98 [8], HEAT [Qj and MASS91 [lO| balloon-borne experiments 
have observed a total of about 80 antiprotons above 5 GeV. However, only two cosmic ray 
antiprotons with a kinetic energy above 30 GeV are reported sl. 

The antiproton-to-proton flux ratio has been measured from 1 to 100 GeV by the 
PAMELA experiment (a Payload for Antimatter Matter Exploration and Light-nuclei As- 
trophysics), a satellite-borne apparatus designed to study charged particles in the cosmic 
radiation with a particular emphasis on antiparticles. The statistics, particularly at high 
energies, is significantly increased compared to the total data sets provided by all previous 
experiments. 

The PAMELA apparatus is inserted inside a pressurized container (2 mm aluminum win- 
dow) attached to the Russian Resurs-DKl satellite and comprises the following subdetectors: 
a time of fiight system (ToF); a magnetic spectrometer; an anticoincidence system (AC); an 
electromagnetic imaging calorimeter; a shower tail catcher scintillator and a neutron detec- 
tor. Technical details about the entire PAMELA instrument and launch preparations can 

n 

be found in jlll]. 

PAMELA has been acquiring data since July ll*'^ 2006. The results presented in this 
letter refer to data acquired in the period July 2006 to February 2008. More than one billion 
triggers have been collected during the total acquisition time of ~ 500 days. Events were 
considered for further analysis if the reconstructed rigidity exceeded the vertical geomagnetic 
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cut-off, estimeted using the satellite position, by a factor of 1.3. Downward-going particles 
were selected using the ToF information. The time-of-fiight resolution of 300 ps ensures that 
no contamination from albedo particles remains in the selected sample. The ionization losses 
(dE/dx) in the ToF scintillators and in the silicon tracker layers were used to select minimum 
ionizing singly charged particles. Furthermore, multiply charged tracks were rejected by 
requiring no spurious signals in the ToF and AC scintillators above the tracking system. 

Particle identification is based on the determination of rigidity by the spectrometer and 
the properties of the energy deposit and interaction topology in the calorimeter. The analysis 
technique was validated using the PAMELA Collaboration's official simulation program 
tuned using particle beam data. 

The tracking information from the spectrometer is crucial for selecting antiprotons. Due 
to the finite spectrometer resolution, corresponding to a maximum detectible rigidity (MDR) 
exceeding 1 TV, high rigidity protons may be assigned the wrong sign of curvature. In 
addition there is a background from protons that scatter in the material of the tracking 
system and mimic the trajectory of negatively- charged particles. In order to accurately 
measure antiprotons, this "spillover" was eliminated by imposing a set of strict selection 
criteria on the quality of the fitted tracks. Track fits required the use of at least 4 (3) 
position measurements along the x (y) direction and an acceptable for the fitted track. 
To remove spillover protons, clean tracking position measurements were required (e.g. no 
accompanying hits from delta-ray emission) and that the MDR, estimated for each event 
during the fitting procedure, should be 10 times larger than the reconstructed rigidity. 
The deflection (1/rigidity) distribution for positively- and negatively-charged down-going 
particles, which did not produce an electromagnetic shower in the calorimeter, is shown 
in Figure [TJ The sample includes events for which the reconstructed MDR is larger than 
850 GV. The good separation between negatively-charge particles and spillover protons is 
evident. As expected, the antiproton tracking requirements limit the distribution of spillover 
protons. 

The calorimeter was used to reject electrons. The longitudinal and transverse segmenta- 
tion of the calorimeter combined with dE/dx measurements from the individual silicon strips 
allow electromagnetic showers to be identified with very high accuracy. Using electrons from 
simulations and particle beams, and simulated antiprotons, we defined an energy dependent 
calorimeter antiproton selection Several topological calorimeter variables are used for 
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deflection (GV ') 

FIG. 1: The deflection reconstructed by the track fitting procedure for negatively- and positively- 
charged down-going particles with a reconstructed MDR > 850 GeV and that did not produce 
an electromagnetic shower in the calorimeter. The shaded histogram corresponds to the selected 
antiprotons. 

the antiproton identification. As an example, the energy density in the shower core weighted 
by the depth in the calorimeter, Qcore/ Nearer is shown in Figure [2l The distribution for the 
proton-dominated positively-charged sample is peaked at 1.25. In the negatively- charged 
sample, the distribution corresponding to electrons peaks at a higher value, and antiproton 
events are collected in a separate peak positioned similarly to that seen in the positively- 
charged sample. The resulting electron contamination was estimated to be negligible across 
the whole energy range of interest. The different, and momentum dependent, interaction 
cross sections for protons and antiprotons were taken into account estimating the calorime- 
ter selection efficiencies as a function of momentum for both species. These efficiencies were 
studied using both simulated antiprotons and protons, and proton samples selected from 
the flight data. In the rigidity interval 2 - 100 GV the proton selection efficiency ranges 
between 0.720 ± 0.003 and 0.800 ± 0.012, whereas the antiproton efficiency ranges between 
0.621 ±0.003 and 0.797±0.012. These efficiencies were used to rescale the number of selected 
antiprotons and protons. 

Possible contamination from pions produced by cosmic-ray interactions with the 
PAMELA payload was studied using both simulated and flight data. Both negatively- 
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FIG. 2: An example of a topological calorimeter variable used for antiproton identification (see 
text for explanation). Positively-charged events are shown in the lower plot. The upper plot shows 
negatively-charged events. The vertical scale for the open histogram has been multiplied by a 
factor of 30 (compared to the filled histogram) for clarity. 

and positively-charged pions below 1 GV were identified using the (3 (velocity) measured by 
the ToF system and the calorimeter information (to reject electrons and positrons). The 
majority of these pion events had hits in the AC scintillators and/or large energy deposits 
in one of the top ToF scintillator clearly indicating that they were the product of cosmic ray 
interactions with the PAMELA structure or pressure vessel. After applying all previously 
described selection criteria, the energy spectrum of the surviving pions was measured below 
1 GV and compared with the corresponding spectrum obtained from simulation by using 



both GHEISHA and FLUKA generators 13|, llJ]. After comparison with the experimental 



pion spectrum below 1 GV, a normalization factor for the simulation, which accounted for all 
uncertainties related to pion production and hadronic interactions, was obtained. The nor- 
malized simulated pion spectrum was used to estimate the contamination in the antiproton 
sample for rigidities greater than 1 GV. This procedure resulted in a residual pion contami- 
nation of less than 5% above 2 GV, decreasing to less than 1% above 5 GV. This result was 
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TABLE I: Summary of proton and antiproton results. 



Rigidity 


Mean Kinetic 




Observed 


Extrapolated 


at 


Energy 




number of 


P 
P 


spectrometer 






events 




GV 


GeV 


p 


p 


at top of payload 


2.23 


- 2.58 


1.64 


39 


1198039 


(3.92 ± 0.63) X 10"^ 


2.58 


- 2.99 


1.99 


48 


1144014 


(4.92 ± 0.71) X 10~^ 


2.99 


- 3.45 


2.41 


55 


1071778 


(5.91 ± 0.80) X 10"^ 


3.45 


- 3.99 


2.89 


60 


988666 


(6.89 ± 0.89) X 10~^ 


3.99 


- 4.62 


3.46 


74 


903708 


(9.2 ± 1.1) X 10"^ 


4.62 


- 5.36 


4.13 


71 


827521 


(9.6 ± 1.1) X 10"^ 


5.36 


- 6.23 


4.91 


93 


738028 


(1.40 ± 0.14) X 10"^ 


6.23 


- 7.27 


5.85 


78 


653736 


(1.31 ± 0.15) X 10"^ 


7 97 


- o.Oo 


6.98 


69 


573172 


(1.32 ± 0.16) X 10"^ 


8.53 


- 10.1 


8.37 


67 


505503 


(1.44 ± 0.18) X 10~^ 


10.1 


- 12.0 


10.1 


94 


449261 


(2.27 ± 0.23) X 10~^ 


12.0 


- 14.6 


12.3 


58 


405583 


(1.54 ± 0.20) X 10"^ 


14.6 


- 18.1 


15.3 


58 


301314 


(2.05 ± 0.27) X 10"^ 


18.1 


- 23.3 


19.5 


46 


270068 


(1.80 ±0.27) X 10"^ 


23.3 


- 31.7 


25.9 


39 


211249 


(1.94 ±0.31) X 10"^ 


31.7 


- 48.5 


37.3 


24 


136858 


(1.82 ±0.37) X 10"^ 


48.5 


- 100.0 


61.2 


6 


57613 


{imt'oiD X 10-' 



cross-checked between 4 and 8 GV by selectin g a ntiproton events below the geomagnetic 
cut-off. This sample includes re-entrant-albedo 23|] antiprotons and locally produced pions. 
By scaling the number of such events for the acquisition time an upper limit for the negative 
pion (and protons with the wrong sign for the reconstructed deflection) contamination in 
the cosmic ray antiproton sample was found to be ~ 3%, in agreement with simulations. 

Table [T] shows the total number of antiprotons and protons that survived the data selec- 
tion. The antiproton-to-proton flux ratio was corrected for the calorimeter selection efficien- 
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cies and for the loss of particles in the instrument itself. It is assumed that all antiprotons 
and protons interacting with the payload material above and inside the tracking system are 
rejected by the selection criteria. The resulting antiproton-to-proton flux ratios are given in 
Table [I] and Figures [3] and HI The reported errors are statistical only. The contamination 
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FIG. 3: The antiproton-to-proton flux ratio obtained in this work compared with theoretical cal- 
culations for a pure secondary production of antiprotons during the propagation of cosmic rays in 



the galaxy. The dashed lines show the upper and lower limits calculated by Simon et al. 



m for 



the standard Leaky Box Model, while the dotted lines show the limits from Donato et al. 16 [for a 
Diffusion model with reacceleration. The solid line shows the calculation by Ptuskin et al. 17| for 
the case of a Plain Diffusion model. The curves were obtained using appropriate solar modulation 
parameters (indicated as 0) for the PAMELA data taking period. 



was not subtracted from the results and should be considered as a systematic uncertainty. 
It is less than a few percent of the signal, which is significantly lower than the statistical 
uncertainty. Figure [3] shows the antiproton-to-proton flux ratio measured by the PAMELA 
experiment compared with theoretical calculations assuming pure secondary production of 
antiprotons during the propagation of cosmic rays in the galaxy. The PAMELA data are in 
excellent agreement with recent data from other experiments, as shown in Figure HI 

We have presented the antiproton-to-proton flux ratio over the most extended energy 
range ever achieved and we have improved the existing statistics at high energies by an 
order of magnitude. The ratio increases smoothly from about 4 x 10~^ at a kinetic energy 
of about 1 GeV and levels off at about 1 x 10""^ for energies above 10 GeV. Our results 
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are sufficiently precise to place tight constraints on parameters relevant for secondary pro- 
duction calculations, e.g.: the normalization and the index of the diffusion coefficient, the 
Alfven speed, and contribution of a hypothetical "fresh" local cosmic ray component 22|. 
Furthermore, an important test criteria for cosmic ray propagation models is their ability 
to reproduce both the antiproton-to-proton flux ratio and the secondary-to-primary nuclei 
ratio. Our high energy data (above 10 GeV) places limits on contributions from exotic 
sources, such as dark matter particle annihilations. The antiproton-to-proton flux ratio will 
be modified according to values of the dark matter particle mass, annihilation cross section, 
and structure in the density profile (boost factor). 

PAMELA is continuously taking data and the mission is planned to continue until at 
least December 2009. The increase in statistics will allow higher energies to be studied. An 
analysis for low energy antiprotons (down to ~100 MeV) is in progress and will be the topic 
of a future publication [13]. 
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(ASI), Deutsches Zentrum fiir Luft- und Raumfahrt (DLR), The Swedish National Space 
Board, Swedish Research Council, The Russian Space Agency (Roscosmos) and The Russian 
Foundation for Basic Research. 
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